Localization of Matter Waves in 2D-Disordered Optical Potentials 
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We consider ultracold atoms in 2D-disordered optical potentials and calculate microscopic quan- 
tities characterizing matter wave quantum transport in the non-interacting regime. We derive the 
diffusion constant as function of all relevant microscopic parameters and show that coherent multiple 
scattering induces significant weak localization effects. In particular, we find that even the strong lo- 
calization regime is accessible with current experimental techniques and calculate the corresponding 
localization length. 
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Introduction: Ultracold atoms in optical potentials 
can be used to realize condensed matter model systems 
in a very versatile manner . Having at hand the pos- 
sibility to shape external potentials almost at will, a 
natural direction of investigation is the disorder driven 
superfluid-insulator transition Q or the Anderson tran- 
sition 0. Interest in this transition has been renewed 
since the experimental observation that even a small dis- 
order in confining fields leads to a fractioning of quasi- 
1D condensates in waveguide structures on atom chips 
. Very recently, several studies of Bose condensates in 
speckle potentials [j| have shown an efficient suppression 
of ID transport by disorder 

The physics of interacting particles in a disordered en- 
vironment has been discussed already for a number of 
years for electrons and superfluid helium 0,0- For cold 
atoms, recent contributions include a numerical study of 
the Bose- Hubbard and Anderson hopping models [10j, 
a renormalization group approach in the case of off- 
diagonal disorder a Bose- Fermi mapping for hard- 
core bosons |l2j . and a transfer matrix treatment of 
atomic matter waves interacting with impurity atoms in 
an optical lattice All these approaches study ID 

systems. Bose condensates in 2D optical quasicrystal lat- 
tices were considered in . In this context, our aim is to 
concentrate on the effect of disorder on quantum trans- 
port in the non-interacting regime, leaving the intrigu- 
ing impact of interactions for later studies. Interaction 
effects are small in the low-density wings of expanding 
Bose condensates Q that are obtained from interaction- 
dominated condensates in the Thomas-Fcrmi regime by 
opening the trap potential. If one insists on working 
in the high-density regime, single-particle dynamics can 
be studied by reducing the scattering length via Fesh- 
bach resonances ^]|, eventually reaching the ideal Bose 
gas regime [lij . Alternatively, one can work with spin- 
polarized fermions 0| whose collisions are blocked by 
the Pauli principle. 

In this Letter, we report analytical results for the 2D 
dynamics of cold atoms in a far-detuned optical speckle 



potential. Using a perturbative Green's function ap- 
proach, we calculate transport quantities relevant for the 
diffusive regime and derive the weak localization correc- 
tion to the classical diffusion constant. The 2D geome- 
try is particularly advantageous for measuring quantum 
corrections to classical transport because diffusive trajec- 
tories always return to their starting point, thus favor- 
ing localization effects. Making use of the microscopic 
characteristics of the speckle potential instead of effec- 
tive models of disorder, we show that a highly disordered 
environment with strong scattering can be tailored with 
speckle potentials. Furthermore, we find that the strong 
localization threshold can be reached with current exper- 
imental techniques. 

Intensity transport: A cloud of non-interacting cold 
atoms is described by the single particle Hamiltonian 
H = p 2 /2m + V(r) where V(r) is a static 2D ran- 
dom potential after the harmonic confinement in the 
transport directions has been switched off as realized 
in [f| 0, 0- The initial state density matrix qq of the 
atomic cloud evolves in time as g(t) = U{t) go U^{t) with 
the evolution operator U(t) = cxp(—iHt/H). Meaning- 
ful matter wave transport observables involve a statis- 
tical average over all possible realizations of disorder. 
One important quantity is the average probability den- 
sity p(r, t) = (r| g(t) |r) of particles at point r and time 
t (a bar denotes the disorder average). Its Fourier trans- 
form p(q, fi) = J d 2 r f dtp(r,t) exp(iilt — iqr) is given 
by 

p(q, Q) = J -0^ Qo (k, q) <&(k, q, Q), (1) 

where g (k,q) = (k + q/2| g |k — q/2) contains all in- 
formation about the initial atomic density distribution, 
and $ is the intensity relaxation kernel for plane waves 
with on-shell energy E = h 2 k 2 /2m. In the long time and 
large distance limits £l,q — ► 0, the 2D relaxation kernel 
for isotropic intensity distributions has the characteristic 
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that describes a diffusion process with diffusion constant 
D(k). Since diffusion soieiy relies on the local conser- 
vation of particles and on linear response it is a very 
robust phenomenon. In the remainder of the paper, we 
will essentially calculate the plane wave diffusion con- 
stant D(k), including interference corrections, as a func- 
tion of all relevant microscopic parameters. Once the dif- 
fusion constant is determined, the dynamics of any initial 
density distribution is obtained by integrating 

Speckle characteristics: Of particular importance are 
the characteristics of the optical potential V(r) — V[l + 
6V(r)] with average V and normalized random com- 
ponent SV(r). The only correlation function that we 
will need in the following is the two-point correlator 
V(r) = <5y(r')<5F(r' +r). We consider two-level atoms 
(transition frequency ujq, transition width T, saturation 
intensity I s ) exposed to a far-detuned monochromatic 
speckle field £(r) at frequency wl = ck^ and detuning 
S = o^l — <-^o- generated from a laser source with power P. 
The speckle pattern is created over a surface of linear size 
L, the local field intensity being I(r) = e$c \£ (r)| 2 /2 with 
average value /l = P/L 2 . The speckle optical dipolar 
potential is V(r) = VI(r)/I L , with V = (hT 2 /8S)I h / I s . 
This random potential derives from the Gaussian ran- 
dom field £ (r), but is not a Gaussian variable by itself. 
Its pair correlator can be expressed as V(r) = |7(r)| 2 
with 7(r) = e c£*(r')£(r' + r)/2iL being the normalized 
two-point field correlation function Its 2D Fourier 

transform "P(k) is the speckle power spectrum. 

For a 2D speckle pattern produced by monochromatic 
illumination of a holographic phase mask (transmission 
geometry) or of a rough surface (reflection geometry), the 
far-field correlation reads 7(1") = 2J\(u)/u where u = r/C, 
and Ji is the first order Bcsscl function [J^- Here 
C = 1/ a^L is the correlation length of the speckle poten- 
tial, and a = R/z <C 1 is the speckle aperture angle at 
a distance z from the speckle source with radius R. The 
correlation length £ defines the intrinsic physical length 
scale of our system. In turn, it also defines an energy 
scale for the atomic dynamics, Eq = h 2 /m£ 2 = 2a 2 Eft in 
terms of the more familiar recoil energy En = h 2 k 2 /2m. 
The ratio r\ = V/Eq measures the strength of the po- 
tential fluctuations relative to the correlation energy Eq. 
The 2D speckle power spectrum is obtained as the con- 
volution of two identical disks, 'P(k) = 8J-(k(/2) with 
T{x) = [arccosa; — xy/l — x 2 \ 6(1 — x). The Heavi- 
side distribution O reflects the fact that the potential is 
smooth on length scales smaller than £, but uncorrelated 
on distances larger than £. 

Weak scattering regime: Microscopic transport pa- 
rameters can be calculated using standard diagrammatic 
Green's function techniques |2(ll2lll22l| . A well controlled 
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FIG. 1: (Color online) Plot of kl s from eq. as a function 
of the reduced matter wave number fc£ for different values 
of the disorder strength n = V/E ( € {0.2, 0.4, 0.8, 1.2, 1.6} 
(thin curves from left to right). The thick line (red online), 
connecting points of kl B where A = 1, i.e. fc£ = V^v, in- 
dicates the limit of validity of the weak scattering condition 
121 The dashed horizontal line corresponds to k£ s = 1/27T. A 
strongly disordered medium with kl B ~ 1 is obtained for cold 
enough atoms. 



perturbative expansion is obtained if the atomic energy 
E lies above the mobility edge E c : 
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In physical terms, this weak scattering condition can be 
understood as a condition for small quantum reflection 
from a potential bump of linear size £ and height V. It 
can be realized in two ways: either the atomic kinetic 
energy E is larger than the correlation energy Eq; then 
the potential fluctuations must be small with respect to 
the atomic energy, V < E. This case corresponds to the 
classical picture of atoms flying well above small potential 
bumps. Or the atomic energy is smaller than the correla- 
tion energy (requiring cooling well below recoil, especially 
if a -C 1); then the fluctuations must still be smaller than 
E^, but can be larger than the atomic energy E. This 
case corresponds to a quantum regime where the atom 
is able to tunnel through high potential bumps of lin- 
ear extension ( thanks to its large dc Broglic wavelength 
A dB = 2ir/k > C- 

Scattering mean free path: In the weak scattering 
limit, the average distance travelled by the atom between 
two scattering events defines the elastic scattering mean 
free path £ s . It is related to the speckle power spectrum 
through 
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where V (k^, 9) = 8J-(k(sin(8/2)) represents the differen- 
tial scattering cross-section. A plot of k£ s as a function 
of the reduced atomic wave number kC, is shown in fig. ^ 
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for different values of the disorder strength r\. The weak 
scattering condition A < 1 implies the bound ki s > l/2n 
such that the shortest achievable scattering mean free 
path is of the order of the 2D speckle correlation length 
£ itself. Fig. ^ shows that the speckle potential, even 
though correlated on the scale f , can become a highly 
disordered scattering medium with kt s of order unity for 
sufficiently cold matter waves (fc£ — > 0). 

Boltzmann diffusion constant: While the atom prop- 
agates through the speckle field, it is scattered by po- 
tential fluctuations. After a large number of scattering 
events, this random walk conserving energy and particle 
number results in diffusive transport for matter waves 
with wave number k. Following the approach pioneered 
by Vollhardt and Wolfle (23, we set up a quantum ki- 
netic equation for the intensity kernel defined in Q that 
allows to calculate the transport mean free path from the 
microscopic properties of the system. In a first step, we 
determine the 2D Boltzmann diffusion constant 
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where the elastic transport mean free path £ tI is the aver- 
age distance travelled by the atom before losing memory 
of its initial direction. This approximation assumes that 
the disorder average washes out all interference effects be- 
tween partial scattered waves. Scattering and transport 
mean free paths are then linked by the relation 



J^dO cos 9 V(k(, 0) 
J^d6V(k(,9) 



(6) 



Because of the potential correlation at small scales (Heav- 
iside function in T(x)), the scattering angle is bounded 
by |sin(0/2)| < l/fc£. Hence, there is no angular re- 
striction for slow atoms (fc£ < 1). In the ultra-cold 
regime (k( -C 1) isotropic scattering ~ is) pre- 
vails, and the speckle potential becomes an effective 8- 
correlated potential. For fast atoms fc£ ^> 1, the max- 
imum scattering angle is max ~ 2/kQ and scattering is 
strongly peaked in the forward direction. In this case, 
£ tr « (kCf) 2 £ s 3> i s (strongly anisotropic scattering). For 
example, for Rubidium atoms with fc£ = 1, L = 2 cm, 
a = 0.1, P = 0.25 W and 8 = 10 6 r, we predict a diffu- 
sive matter wave transport with elastic scattering mean 
free path £ B ks 2 /im and elastic transport mean free path 
£tr ~ 7 /im. 

Weak localization: In phase coherent samples, the 
constructive interference between counter-propagating 
amplitudes enhances the return probability of the atomic 
matter wave to a given point. This weak localization 
correction (22l (^ reduces the diffusion constant to D = 
D B - SD with 
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FIG. 2: (Color online) Weak localization corrections 8D rel- 
ative to the Boltzmann diffusion constant Db at detuning 
8 = 10 6 r, eq. 0, as a function of laser power P for differ- 
ent atomic matter wave numbers (from left to right) fc£ = 
1.25, 1.5, 1.75, 2.0, 2.25, 2.5. The speckle size is L = 2 cm, 
the aperture angle a = 0.1. For each value of the weak 
scattering condition A < 1 is valid to the left of the thick 
black line (solid curves). Already for moderate laser power, 
weak localization corrections of observable size are predicted. 



in 2D. Here, the length Lq = min(L, L^) is the relevant 
cutoff for fully coherent multiple scattering. It is cither 
the system size L itself, or the phase coherence length 
L<t> = \/DBTi which accounts for possible phase break- 
ing mechanisms affecting interference at a rate Ti = rf 1 . 
For cold atoms in optical speckle potentials, one phase 
breaking mechanism is inelastic scattering of photons for 
which Ti = TV/H8. The corresponding phase coherence 
length Lfi scales as 8 2 I^ 3 ^ 2 while i s and i tr scale as 8 2 l£ 2 . 
Thus, by keeping the ratio Il/8 fixed, all multiple scat- 
tering parameters V, £ s , itr and Db are kept constant. 
Under this condition, changing 7l (and 8) only modifies 
Ltj, and thus the interference corrections. This opens the 
way to use inelastic scattering to monitor the weak local- 
ization corrections in a controlled manner, just like with 
an external magnetic field in 2D electronic experiments 
on negative magnetoresistance |24j . 

In fig. |21 we plot 8D/Db as a function of the total laser 
power P for different values of the atomic matter wave 
number k at fixed laser detuning 8 = 10 6 T. The thick 
black line indicates the corresponding limit of validity of 
the weak scattering condition A < 1. The colder the 
atoms, the larger the quantum corrections (within the 
lower bound £/L for fc£ and for the applied laser power 
imposed by the diffusion condition Lq > £ ir ). Even mod- 
erate laser power assures coherent multiple scattering in 
the speckle plane and induces sizeable weak localization 
corrections. For fc£ = 2, the relative correction SD/Db 
attains approximately 30% at P = 0.9 W. For kQ = 1.25, 
the border 8D = Db is even reached within the region 
A < 1 at P = 0.5 W: this is the strong localization onset. 
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to the strong localization regime is experimentally ac- 
cessible and therefore worthwhile further numerical and 
experimental investigation. We have given an estimate 
for the localization length and provided a set of reason- 
able experimental parameters that we hope will facilitate 
the realization in the laboratory. 
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FIG. 3: (Color online) Logarithmic plot of the 2D localization 
length £io C and of the phase coherence length as a function 
of the laser power P for fc£ = 1.25 and 8 — 10 6 Y. The system 
size is fixed at L — 2 cm, the aperture angle at a = 0.1. £i oc 
and L$ cross at the strong localization threshold, which is 
reached for P — 0.5 W, where the corrected transport mean 
free path £™ L = 2m(DB — SD)/hk vanishes. 

Towards strong localization: The perturbative weak 
localization correction Q diverges with the cutoff length 
Lq, which is compatible with the scaling prediction that 
waves are always localized in 2D on sufficiently large co- 
herent length scales [25| . Equation J7J) is actually the re- 
sult of a self-consistent theory [2j| that takes into account 
the relevant singular terms driving the system towards 
the strong localization onset SD/D-q — > 1. The strong 
localization threshold is then reached when Lq reaches 
the localization length £i oc ~ t s cxp(7rfc^ tl /2). 

Because of its exponential growth, the localization 
length becomes very large when k£ s increases. Observing 
strong localization requires fully coherent scattering in a 
sufficiently large speckle field. We therefore have to check 
whether these requirements can be met with reasonable 
laser power and atom temperatures. In fig. [21 we have 
plotted £io C and as a function of P for fc£ = 1.25, 
L = 2 cm, a = 0.1, and <5 = 10 6 r. The two curves cross 
when P = 0.5 W. At this point, we find for Rb 87 atoms 
(Al = 27r/fcL = 0.78 /xm) £i oc = 2 mm, £t t = 5 /xm and 
£ s = 0.92 /xm. This places the strong localization thresh- 
old at kt s ~ 0.93 and -q sa 0.77, meaning that the atoms 
have an energy slightly above the speckle fluctuations, 
E fn 1.02 V. The atomic temperature is then T = 2.8 nK 
which is experimentally accessible |2fij . The correspond- 
ing de Broglie wavelength is AdB = 8A^ sa 6 /xm. 

To summarize: Making use of quantum transport 
theory, we have studied the dynamics of ultracold atoms 
in 2D speckle potentials. Starting from the microscopic 
potential correlation function, we have calculated the 
elastic scattering mean free path and the classical dif- 
fusion constant. Weak localization corrections are shown 
to be of measurable size for realistic laser power and 
atom temperature. We have found that the threshold 
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